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Abstract 

The localization technique allows us to compute the free energy of the U{N)i. x 
C/(A^)„fc Chern-Simons-matter theory dual to type IIA strings on AdS4^ x CP^ from 
weak to strong 't Hooft coupling A = N/k at finite N, as demonstrated by Drukker, 
Marino, and Putrov. In this note we study further the free energy at large 't Hooft 
coupling with the aim of testing AdS/CFT at the quantum gravity level and, in particu- 
lar, sum up all the corrections, apart from the worldsheet instanton contributions. 
The all genus partition function takes a remarkably simple form - the Airy function, 
Ai((7rA;^/\/2)^/^Aren) ) with the renormalized 't Hooft coupling Aren- 
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1 Introduction and main results 



There has been considerable progress in testing the AdS/CFT conjecture [T] in the large 
limit. Through its connection to the spin chain systems, the large integrability allowed us 
to successfully check the duality from weak to strong 't Hooft coupling with remarkably high 
accuracy [2]. In the meantime, there is hardly any test of the AdS/CFT duality at finite A^. 
On the one hand, the 1/A^ corrections are quantum gravity loop corrections and thus very 
hard to compute on the gravity side. On the other hand, there are indications from the field 
theory analysis that the integrability breaks down at the non-planar level [5]. So there was 
virtually no technique to evaluate the 1/A^ corrections at strong 't Hooft coupling. 

However, there has been important progress in the localization technique in supersym- 
metric gauge theories [1] , elucidating ealier works |5l |6] . This provides a powerful method 
to compute the 1/A^ corrections at strong 't Hooft coupling in dual CFT's. In [7j the lo- 
calization technique was applied to A/" = 6 U{N)k x U{N)-k Chern-Simons-matter theory, 
also known as the ABJM theory, dual to type HA strings on AdS4 x CP^ or M-theory on 
AdS^ X S"^ /Zk [8j. The partition function of the ABJM theory reduces to the eigenvalue inte- 
grals and this defines the ABJM matrix model. Subsequently, the free energy of the ABJM 
matrix model was computed exactly in the large A^ limit and order by order in the 1/A^ 
expansion P [10, [H]- In particular, the genus zero free energy at strong 't Hooft coupling 
precisely agrees with the tree level SUGRA action. This provided a remarkable check of the 
AdS / CFT duality and showed the power of the localization techniquejl] 

In this note we study further the free energy of the ABJM theory at strong 't Hooft 
coupling at finite A^ with the aim of testing AdS/CFT at the quantum gravity level and, in 
particular, resum all the 1/A^ corrections. Our main result is that, apart from the worldsheet 
instanton contributions, the all genus free energy of the ABJM matrix model sums up to 
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where Ci = {lix jk)-^!'^ j ^ and A = Njk is the 't Hooft coupling of AT = 6 V{N\ x f/(A^)_fc 
Chern-Simons-matter theory, and 

1 

Arcn = A - ^ - ^ • (1.2) 

This non-planar shift / renormalization of the 't Hooft coupling (11. 2p partially resums the all 
genus free energy [TT]. As we will see, once this partial resummation is done, the remaining 
higher genus free energies obey a very simple recursion relation [16] which one can easily 
solve to find the main result (II. ip . 



^See related interesting developments; generalizations to quiver gauge theories |12| . the i^-maximization 
conjecture |13 , computations of superconformal index [Tl], and tests of three-dimensional dualities [T5] . 
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It is worthwhile to note that the shift (11 ■2p is closely related to the renormalization of 
the AdS radius (in the a' = 1 unit) [T7] : 



In the planar limit N oo, the matrix model shift (II. 2p agrees with that of the dual string 
theory (II. 3p . as noted in However, there appears to be a discrepancy at the non- 

planar level; The matrix model suggests a further shift — 3A^/8A^^ to the string/M theory 
prediction. This may provide a simplest test of AdS/CFT at the quantum gravity level, and 
it is very important to understand whether and how this discrepancy can be reconciled. We 
will return to this point later in the discussions. 

The organization of the rest of our paper is as follows: In section H] we briefly review the 
localization technique and the ABJM matrix model and summarize the result of [H] relevant 
to our study. In section [3] we give a brief review of the technique to compute the higher genus 
free energies by using the so-called holomorphic anomaly equation and carry out a partial 
resummation of the all genus free energy. In section H] we propose the all genus free energy 
up to the worldsheet instanton corrections and give a proof of our proposal. In section we 
summarize our results and discuss the SUGRA one-loop corrections to the free energy. In 
the appendices we give the technical details of our computations. 



The localization technique [1] allows us to compute the free energy of the ABJM theory from 
weak to strong 't Hooft coupling A = N/ k at finite N. The localization of this theory on 
was carried out by Kapustin, Willet, and Yaakov in [7]; The partition function reduces to 
the finite dimensional integral over the eigenvalues of two U{N) matrices. This defines the 
ABJM matrix model. As demonstrated by Drukker, Marino, and Putrov [9l[TTl[T0], the free 
energy can be calculated for arbitrary 't Hooft coupling A and order by order in the 1/N 
expansion. 

The field content of the ABJM theory consists of two 3d A/" = 2 U{N) vector multiplets, 
{A^, a^, A^, A"^, D^)a=i,2, which are the dimensional reduction of 4d A/" = 1 vector multiplets 
and two bifundamental chiral multiplets, (0^, 0^, V^^, '?/^'^, F^, 2, in the representation 

{N,N), and their duals in (iV, A^)!^] This theory enjoys the nilpotent Grassmann-odd sym- 
metry generated hj 6 = eQ where e being the two-component complex Killing spinor on 5*^ 
with the normalization ee = 1 and Q is the Af = 1 supercharge. 

^In [5 the unequal rank U{Ni) x U{N2) case [18] was also studied, and the generalization of (jl.2|) precisely 
agrees with the SUGRA prediction [19] in the large N limit. 

■^In Euclidean space, the barred fields are not complex conjugate of the unbarred fields but independent 
of them. 




(1.3) 



2 The ABJM matrix model 
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The partition function is invariant under the deformation of the action by any 5-exact 
terms, since the A/" = 6 U{N)k x U{N)_k Chern-Simons-matter action and the vacua are 
invariant under the 5-transformation. This has an important ramification: By deforming the 
action by t6V with a positive definite 6V, as we send t to infinity, the path integral locahzes 
on the saddle point of 6V. A particularly convenient choice of 6V = SVgange + ^Knatter is 

= ^S^^uiN) (^^^^^ - = ^YM , (2.1) 

^Kxattcr = M {^/j'tfj' - 2t^' {a' - a^) <j)' - = 5m , (2.2) 

where Sym and are the N = 2 super Yang-Mills action and the bifundamental matter 
action on the unit sphere, respectively. Their saddle points are given by [7] 

A;^ = 0^ = , = -o^ = const . (2.3) 

The localization action S'ym + 5'm vanishes on the saddle points. Thus, in the t — i- oo 
limit, the only contribution comes from (1) the classical Chern-Simons-matter action Sd 
evaluated on the saddle points and (2) the quadratic fluctuations about the saddle points in 
the localization action. The latter gives the one-loop determinants, and thus schematically 
the partition function becomes 

Z-"^/<;.'d.=|i|£gexp(-S.,(.-^)). (2.4) 

where and Ap denote the Laplacians on S'^ for the bosonic and fermionic fluctuations, 
respectively. 

Diagonalizing the U{N) matrices a*^^'^ and integrating their angular parts yields the 
Vandermonde determinants which are cancelled by the factors from the one-loop determi- 
nants. The net result is given by [3 [221 El [ID] 

.ABJM. ^ / j^^j^ci..nK.(2smh(^))^n.<.(2sinh(i^))^ 



(^0^ J t\ 2^ t\ 2- n., (2 cosh {l^)Y 

X e-5i7(^»'^?-S-^') , (2.5) 

where Qs = 2Tii/k. /Xj's and i/^'s are the eigenvalues of the two U{N) matrices (Xj, and the 
hyperbolic functions are the one-loop determinant contributions. Those in the numerator 
are from the vector multiplets, whereas those in the denominator are from the bifundamental 
matter multiplets which give the coupling between the two U{N) factors. The eigenvalue 
integral (12.51) defines the ABJM matrix model. 

The planar limit of the ABJM matrix model can be solved by the standard technique 
[SI [m [ID]. In the large N limit, the partition function (12. 5 p is dominated by the saddle 
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points: 

/i, = ^ y coth ^^^^ + tanh ^^i^ , (2.6) 

. N2 . 

fe^^a ^ i=l 

where ti = QgNi and A^^i = —N2 = N. In solving these equations, one first considers the 
problem for positive A^^i and N2 and then analytically continues the result to the negative 
N2 = -Ni = -N. 

The effective potential which leads to the above saddle point equations consists of attrac- 
tive harmonic potentials and logarithmic Coulomb repulsions for the eigenvalues /Xj's and 
Ua + ivr's. Since the Coulomb repulsions have the strength tj's, when the 't Hooft couplings 
tj's are zero, the harmonic potentials dominate and the eigenvalues collapses to zero. As we 
increase tj's, the Coulomb repulsions kick in and spread the eigenvalues along two line inter- 
vals; The /ij's condense in Ci on the real axis, and the z/^ + ivr's in C2, the interval separated 
from Ci by +iTT. Thus the ABJM matrix model is a two-cut model. 

Similarly to the case of the standard matrix models, one defines the resolvent by 

c^W = ^7s/5^^°*^^Y^)+^-(5^*^^h^T^) • (2.8) 

To make more direct contact with the standard matrix models, it is useful to introduce the 
new variable Z = e^. Then the resolvent is expressed as 

where t = ti + 12. As studied in [21], the saddle point equations (12. 6p and (12. 7p imply that 
the discontinuity of the resolvent ujo{z) = limAr_j.oo i^(-2) at large is given by 

z = - {uJo{z + ie) + ujo{z ~ ie)) on Ci , (2-10) 

z = - {uJo{z + iiT + ie) + Uq{z + iiT — ie)) on €2- (2-11) 

It then follows that the function 

f{Z) = e* (e"^" + Z^e"^") (2.12) 

is regular everywhere on the complex plane. This has the asymptotic behavior f{Z) ^liS? 
and f{Z) — > 1. With these boundary conditions, the function f{Z) is uniquely determined: 

f{Z) = Z'-CZ + l . (2.13) 
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Meanwhile, from fl2.12p one has 

Uo{Z) = log 



(2.14) 



Parameterizing the branch cuts Ci and C2 on the Z-plane by [1/a, a] and [—6, —1/6], respec- 
tively, the inside of the square root can be written as the polynomical {Z — a){Z — l/a){Z + 
b){Z + 1/6), and one can then identify 

Once the resolvent is found, the 't Hooft couplings can be computed from 

= ujoiz)dz , (2.16) 

and the genus zero free energy from 

= Tcit - ^ ^ Uo{z)dz , (2.17) 

where s = |(ti — ^2)- The contour Ai encircles the branch cut Ci, whereas the contour B is 
the cycle dual to the cycles. The first term in the RHS of (12.1 7p is due to the pole at 
Z = in (12. 9p . These contour integrals (I2.16P and (I2.17P are hard to carry out. However, 
once the derivative is taken w.r.t. (, it becomes easy to perform the integration. 

After the analytic continuation ^2 — > ~ti = ~27riy = 27rzA, the answer for the ABJM 
matrix model turns out to be [HI [TH 
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where k = —iC,. In particular, at strong coupling A 1, these yield 

A - 1 = 1^ + 0(l/«:2) , (2.20) 



24 27r2 

F„(A) = M^! (a - i) O (e--V^) . (2.21) 



Quite remarkably, the first equation shows that, at least in the large N limit, the free energy 
of the ABJM theory at strong coupling is naturally expressed in terms of the shifted 't Hooft 
coupling A — ^, rather than A, in agreement with the SUGRA prediction [17]. Moreover, 
the first term in the genus zero free energy precisely agrees with the (minus of) classical 
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SUGRA action. As A = N/k^ this in particular reproduces the tantahzing N^/"^ scahng of 
the M2-brane theory. Lastly, the exponential corrections to the free energy agree with the 
expected worldsheet instanton corrections. 

These results proclaim the power and relevance of the localization technique to the test of 
the AdS/CFT conjecture. With this remarkable success, it is natural to ask if the AdS/CFT 
can be tested beyond large limit by studying the non-planar corrections in the ABJM 
matrix model. 

3 Higher genus free energies 

The ABJM matrix model is equivalent to the so-called Lens space matrix model [231 121], 
which computes the partition function of the U{Ni) x U{N2) branch of the Chern-Simons 
theory on the Lens space L(2, 1) = S'^ /Z2-, by the analytic continuation A'2 — )■ —N2 [2^ . 

Meanwhile, as discussed in plj, the Lens space matrix model arises as the low energy 
effective theory of D-branes wrapping S^/Z2 in the topological A-model on T*{S^/Z2). At 
large the geometric transition takes place and S^/Z2 is replaced by the Hirzebruch surface 
Fq = X P^, yielding the topological A-model on a non-compact Calabi-Yau space, the 
canonical line bundle over Fq. This is the large A^ (open/closed string) duality between 
the Lens space matrix model and the topological A model on local Fq. Furthermore, the 
A-model on local Fq can be mapped to the B-model on the mirror manifold given by the 
surface uv = H{x,y), where H{x,y) = is an elliptic curve and coincides with the spectral 
curve of the Lens space matrix model. The closed string amplitudes F^^\t,i) of the B- 
model, where t is the complex structure moduli, obey the holomorphic anomaly equation 
of [25] and can be explicitly calculated. In [26] it was conjectured that the holomorphic 
limit F^9)(^t'j = \imt^^F'^3\t,i) of the amplitudes are given by the genus g free energies of 
the matrix model whose spectral curve is H{x, y) = 0. Direct proof was given in [27] built 
on [28] that the solution of the matrix model loop equation [29j is that of the holomorphic 
anomaly equation in the t — t- 00 limit when the appropriate boundary conditions are imposed. 
Thus one can use the holomorphic anomaly equation to find the genus g free energies of the 
ABJM matrix model which is equivalent to the Lens space matrix model by a simple analytic 
continuation. 

3.1 The holomorphic anomaly equation 

Parameterizing the complex structure moduli by the coordinates , the holomorphic anomaly 
equation takes the form [23]: {g > 2) 

dfFg = hojjKe^''G''G''^ (^jDkF^_, + j , (3.1) 
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where Gjj = djdjK is the Kiihler metric on the moduh space, Cjjk is the Yukawa couphng 
defined by Cjjk = djdjdxFQ with Fq being the genus free energy, and the covariant 
derivative Dj acts on Fg and its derivatives as dj — F*, + (2 — 2g)djK. In the local Calabi- 
Yau case, the Kahler potential is given hj K = ^ (t^dj^pQ — t^OxFo) , and thus the moduh 
space metric is Gjj = Imr/j with tjj = djdjFo 



When applied to the ABJM matrix model [9], the moduli space becomes one-dimensional 
and the modulus t is identified with the 't Hooft coupling A. Then the Yukawa coupling is 
Cxxx = 4a|Fo(A) = -327r3z^ where 

with the normalizations properly adjusted. The genus g free energies are assumed to be of 
the form 

f,{t) = e'-'m , (3.3) 

Ur^-E^mU^irU (3.4, 

k=l ^ ^ 

where Af^\T)^s are the modular forms of weight 6g — 6 — 2k, and c^(r) is the holomorphic 
ambiguity and a modular form of weight 6g — 6. The Eisenstein series of weight 2, E2{t), is 
a quasi-modular form and can be promoted to the non-holomorphic modular form as 

E^ir) ^ E^ir, f ) = E,{t) - . (3.5) 

Imr 

By promoting Fg{T) to the non-holomorphic modular form Fg{T,f) in this way, the non- 
holomorphicity comes in only through £"2. It then follows from the holomorphic anomaly 
iKT} that {g > 2) 

^ = l^/.-! + l^d^f^-^ + E Wc/.-.} , (3.6) 

where the covariant derivative transforms a form of weight kto a form of weight k+2: d^ = 
9,- + In Appendix A we collect various formulae used in the subsequent computations. 

The genus one free energy can be independently calculated by using Akemann's formula 
[30] . One can find that 

Fi (r) = - log r7(r) , (3.7) 

where 77 (r) is the Dedekind eta function. This is the initial data needed to solve the holomor- 
phic anomaly equation. In particular, at strong coupling A ~ log^ k ^ 1, this is expanded 

1 1 

Fi(r) = - log - - log (log k) + 0{1/k) . (3.8) 

D 2 
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3.2 The weight zero free energy 

We first focus on y4„ := A^g^ ^■''s of modular weiglit zero. By using tlie formulae in Appendix 



A, it is easy to show that they obey a simple recursion relation [16 

A 



'■9 

12 



^(3^7-3)A,_i + ^A,Vr j • (3.9) 



They are the highest order terms in F^'s in the x = 1/log/t expansion. Indeed, using the 
modular transformation properties in Appendix A, one finds that 

Fg = + O {x'^-') := Ff x^^"^ + O {x'^-') , (3.10) 

where x = 1/log/t. As it will turn out, the leading order recursion (13. 9p carries almost all 
the information needed to sum up the all genus free energy, when the worldsheet instanton 
corrections are ignored. 

We can in fact resum the highest order all genus free energy 

oo 
3=2 

where t = —iggX^^"^. For the convenience, we instead consider the generating function 



oo 



/(«):= (3.12) 

9=2 

where Ag = {g — 1)(4/3)^~^fJ°^. Then it is easy to find that 

f{u)-A,u' = u'f{uy + f{uf (3.13) 
with A2 = 5/36 which can be found from the initial genus one data. This can be solved to 

fiu) = ~^\ogHiu) , (3.14) 

where 

Hiu) = e-^y^T/2^{CiKi/^s{l/2u) + C2li/3{l/2u)) , (3.15) 

with Ci and C2 being integration constants. The highest order all genus free energy can 
then be computed as 

F{t) = - du-^f{u) = logi/(3tV4) • (3.16) 

Note that we absorbed the integration constant into Ci and C2. 
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To fix the constants Ci and C2, we require the absence of the non-perturbative corrections 
of the type e"^/^*^ for small t. Since t ^ Qs, this is of the order e~^/^» and would be the 
gravitational instanton or NS5-brane effect, as opposed to the D-brane effect, in the dual 
AdS string theory. We assume that the type IIA string theory on AdS^ x CP^ does not 
receive such corrections. With this assumption, there remain two choices; (1) C2 = or 
(2) Ci = — (— l)^/^C2/7r. However, the latter choice is subject to a further non-perturbative 
ambiguity due to the Stokes phenomenon of the modified Bessel function Ii/^^z). So we 
claim that the correct choice is C2 = 0. The remaining constant Ci only shifts the constant 
part of the genus one free energy and is not physically important. 



Using the identity Ai(z) 



i^i/3(i^=^/'),wefind 



F (t) = log 27iCie~^t~3Ai{t' 



-4/3N 



(3.17) 



where t = —igsX^^'^. Note that the factors in front of the Airy function are precisely the 
minus of the highest order terms (in x) of the genus zero and one contributions. Hence the 
total weight zero free energy becomes simply 



{gs,x) = \og 



27rC,Ai({-gy) 



-2/3 



(3.18) 



As we will see momentarily, we have actually resummed a significant part of the all genus 
free energy up to the worldsheet instanton corrections. 



4 The all genus free energy 



It was observed in [TTj that the all genus free energy, apart from the worldsheet instanton 
corrections, depends on the 't Hooft coupling only through the shifted/renormalized form 

This translates to the renormalization of QgX by 

OgX — )• = — , (4.2) 

where Qs = 2mX/N and = v^tt^A - 1/24. 

We now claim that, apart from the worldsheet instanton corrections, the all genus free 
energy can be obtained by replacing x with the renormalized variable y in the weight zero 
free energy fl3.18p . That is, the all genus free energy of the ABJM theory is given by 

^ABJM ^) ^ j27rCiAi (^{-gyy'^'^j + O , (4.3) 
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where Ci = {—g1) /\/2 to match the normahzation in yULj. This is the main result f ll.ip . 

To show it, we first note that, when the worldsheet instanton corrections, i.e., the 0{1/ k) 
terms, are neglected, the holomorphic anomaly equation (13 .Op becomes 



I2x - 1 



-X 



9-2 



X) 



(4.4) 



r=2 



To derive this, we performed the modular transformation ST^^ : r i— ?► r' = — l/(r — 1), and 
kept only the terms in powers of log/t = -nir' /2 + O (e^'^*'^'). We then used the formulae 
listed in Appendix A. 



To proceed, we define the generating function 



g=2 



X) 



(4.5) 



which obeys 



1 



12a; - 1 



J^{9s, x) - glF!^{x) = -glx^d^T{gs, x) + 



gsX^J^igs, x) + ^g^,x^J^{gs, xf 



where 



h^ix) = \ 

^ ' 144 24 48 



(4.6) 



(4.7) 



Then the claim (14. 3 p is equivalent to 



J^{gs,x) 



dy / 4 

dx \ gly"^ 



f ( ] - d, 



n 1, 

+ - — h - log X 



where 



3ggX^ \6x 2 

1 z' 
-2 + 6 



(4i 



(4.9) 



is the generating function (13.141) plus the genus zero and one contributions, while the second 
term in (14. 8 p is the subtraction of the genus zero and one contributions. Plugging (14. 8 p into 
(iSD yields 



/ ( 



36 V 4 



gy ) /' ( -yy ) + / ( -yy 



(4.10) 



This is equivalent to (I3.13p . Hence we have shown that the all genus free energy (14.31) indeed 
satisfies the holomorphic anomaly equation (13.61) up to the worldsheet instanton corrections. 



11 



5 Discussions and conclusions 



We succeeded to sum up the all genus free energy of the ABJM theory, apart from the 
worldsheet instanton corrections. The resummation was done in two steps; (1) First, we 
summed up the modular weight zero part of the all genus free energy as given in f l3.18p . (2) 
Then the remaining parts were resummed to (14. 3 P by renormalizing the 't Hooft coupling as 
was done in [TT]. The renormalization (11. 2p agrees with the SUGRA prediction [17] in the 
large N limit. However, there seems to be a discrepancy at the non-planar level; The matrix 
model suggests a further shift — 3A^/8iV^ to the string/M-theory prediction (II. 3p . 

The quantum gravity one- loop shift +A^/24A^^ of the AdS radius in (II. 3p accounts for, at 
the least, a part of the one-loop corrections. This comes from the shift of the M2/D2-brane 
charge [17] 

N^N- — (k-l-^=k(x- — + , (5.1) 

24 V kj V 24 24iV2 y ' v ; 

due to the higher curvature correction C3 A Jg, where /§ is a curvature 8-form anomaly 
polynomial [3T]. One might wonder if the one-loop discrepancy — 3A^/8A^^ comes from other 
SUGRA one-loop corrections such as the R'^ term [32] E However, we now argue that they 
are absent: In fact, the one- loop R'^ term is nonzero and of the order a/A which agrees with 
the leading part of the genus one free energy in (13. 8p . However, the numerical factor does not 
agree with the — 3A^/8A^^ discrepancy, and the i?^ term is not the only one-loop correction. 
The one- loop i?^ term is expected to be completed by the terms involving the 4- form F4. 
Indeed, the lid J\f = 1 SUGRA on-shell superfield [33] suggests schematically the completion 
of the form {R + -Fl)"^, but this vanishes on AdS^ x S'^ /Zk [SUES]. Furthermore, the one-loop 
vacuum energy vanishes in 4d A/" > 5 SUGRA on the global AdS^ whose boundary is M x S*^ 
|36j . Although we work in the Euclidean AdS^ with the boundary S^, the localizations of the 
ABJM theory on and Mx 5*^ turn out to be the same [37] . So we expect the SUGRA result 
on the global AdS^ applies to our case, and it seems likely that the one-loop corrections, 
except for the charge shift (15.11) . are absent. 

This leaves us the — 3A^/8A^^ discrepancy and suggests us to look for the resolution on 
the matrix model side. Recall that ABJM proposed two A/" = 6 superconformal Chern- 
Simons-Matter theories [8]; One has the gauge group U{N) x U{N) and the other has 
SU{N) X SU{N). In the large N limit, their difference may not matter. But at finite N 
they will differ from each other in the 1/A^ corrections. In fact, in the AdS^ x 5*^ case, the 
correct gauge group is SU{N) rather than U{N) [38]. Although the status of this subtlety 
is unclear in the AdS4 case, it may be worthwhile to investigate SU{N) x SU{N) ABJM 
matrix model to see if the — 3A^/8A^^ discrepancy can be resolved. We hope to address this 
issue further in the near future. 

''The \^/N^ shift appears of the order \/A in the genus one free energy. 
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Finally, it is worth emphasizing the remarkable simplicity of our result; The partition 
function of the ABJM theory, when the worldsheet instantons are neglected, is simply the 
Airy function. This might suggest a possible connection of the ABJM theory to the Kont- 
sevich matrix model [39] upon the inclusion of the worldsheet instanton corrections. Mean- 
while, in a somewhat different context of M-theory flux compactification, it was argued that 
the norm square of the "wave-function of the universe" or the 5d black hole entropy was 
given by the Airy function [40j It would be interesting to study the relations to these works 
in the future. 
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A Formulae for holomorphic anomaly equation 

In this appendix, we gather various useful formulae related to the holomorphic anomaly 
equation: 



,4 



^ " ' ^^•^^ " ' d^b=^b{b + d), d(.d = -^d{b + d) , 

The following formulae are useful to derive the weight zero recursion (13. 9p 

^(35-3) / \ 3g-3 



3(7-3 vi2y 



if) +■■■• ^^-2) 



To study the perturbative terms as in (14. 4p . we shall use the modular transformations ST ^ : 
r I— 7- r' = — 1/(t — 1), since (13. 2p implies that r is given by r = 1 + 7rz/(21ogK), neglecting 



^The authors are grateful to Hirosi Ooguri for pointing this out to us. See also 
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the 0{1/k,) instanton terms, 



biT) = -{-^Trd{r') = {Tr[l + Oil/^)] , 

d{T) = {-ZT'yc{r') = {r')'[-l + 0{l/^)] , 
E,{t) = (r')^[i?2(r') - Q^/{nT')] = {r')\l - 3x) + (9(1//.)] , 
e{T)=4/{b{r)d{r)') = {rr'[^ + 0{l/K)], 



(A.3) 



with X = 1/ log K. Note that since the holomorphic anomaly equation is a modular covariant 
equation, the modular scale factors r' are destined to be cancelled out in the end. We 
shall drop them in our following analysis. Then, we obtain the following formulae in the 
approximation neglecting the 0{1/k) terms: 



d^E2 



1 dr^ X 



dE2 
dx 



-3 , d^b = d^d = . (A.4) 



From these we find 

dh 



dEo 



difg = {diE2, 



dfg 

dE2 

dfg 



X 



4 ^^(-) 



difg = {diE2)^ + (d^E, 



dE^ 



S 4 



(A.5) 



B A derivation of weight zero free energy 



In this appendix we solve the differential equation flS.lSp : Introducing the new variable 
V = l/{2u) and the new function 



this becomes 



d^ A2 

h 2 \ 

dv"^ dv 



H = . 



(B.l) 



fB.2) 



In terms of G{v) = e'"v ^H(v), one obtains (a = 1/3) 



2d d / 2 I 2\ 



G = 



(B.3) 



which is nothing but the modified Bessel's differential equation. The solution is given by the 
modified Bessel's functions: G{v) = CiKi/^^v) + C2h/3{v). 
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